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Abstract 

We discuss the implications of the recent evaluations of the SM contribution to a^ = (g- 
2)^/2 in the light of the latest E821 measurement, which indicate ~ 3a deviation. We 
derive the 95% CL interval, Sa^, to be used to constrain any additional contribution to 
a^, beyond the SM ones; it has to have a positive sign. We apply the 5a ^ to constrain 
the light Higgs-boson scenarios in a Two-Higgs-Doublet-Model ("Model II"). If the con- 
straints from the new (g — 2)^ results are combined with the other existing constraints, 
one can exclude a light-scalar scenario at 95% CL, while a light-pseudoscalar scenario 
can be realized, for a pseudoscalar mass between 25 and 70 GeV, with tan j3 in the range 
25 < tan/3< 115. 
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1 Introduction 



A precision measurement of the g — 2 for the muon at BNL is expected to test the electroweak 
(EW) sector of the Standard Model (SM) and at the same time to shed light on possible effects 
of "new physics". After a release of the new E821 result |], based on the n + data collected in 
the year 2000, a current mean of experimental results for (g — 2)^ is (Q: 



( a — 2) ex P 

a ex P = g ! ± = n 65g ^ x ^-10^ ^ 

with an uncertainty (in parentheses) that is almost two times smaller than in the previous 
measurement ||, and only two times larger than the ultimate goal of the E821 experiment. 

The Standard Model prediction for consists of the QED, EW and hadronic contributions: 

= a Q * D + a™ + a™. 

F P p P 

The QED contribution, which constitutes the bulk of the SM contributions, is calculated up to 
four loops and the 0(a 5 ) term is estimated ||; its uncertainty is very small ~ 3 x 10~ n . The 
EW contribution, based on one- and two-loop diagrams, is also known to a similar accuracy 
[4-8] The EW contribution is small (152 x 10~ n ) and is only about two times larger than 
the present experimental uncertainty (see (1)). The hadronic contribution, ~ 7000 x 10~ n , 
is the second largest contribution to a^ M . It is known at present with an accuracy of 1 %. 
This uncertainty is the main source of the uncertainty in the SM prediction. Predictions for 
the hadronic contribution [9-21] differ among themselves, however, these differences seem to 
be much less significant than, say, one year ago (see e.g. discussions in [11,12,17-21], and also 
in [E3L El]). The dominant contribution to a^ ad , as well as one of the dominant errors in its 



value, come from the leading vacuum-polarization (vpl) term. Some preliminary results of the 
improved calculations of the vpl have been presented recently []T3], [TJJ. These are data-driven 



analyses using the most recent data on hadron production in e + e~ collisions from BES, CMD-2, 
SND |p2| . The corresponding uncertainties for the vpl, ~ 58 x 10~ n , are now even smaller 
than those obtained previously by using the data on r decays in addition to the then available 
e + e~ data, e.g. see (l^, [l^]. Another important issue has been the hadronic contribution to the 
light-by-light (lbl) scattering, and its contribution to (g — 2)^. After a sign error, first pointed 
out in was found in the earlier calculations of this contribution, a few re-evaluations Jl9| 
of this part have appeared during the last few months. All of them confirm the finding of 
[ |j~8|j , i.e. a positive sign of the lbl contribution. The central value for the lbl varies from 80 to 
110 x 10~ n , depending on the analysis. Since this contribution can be estimated only on a 
purely theoretical ground, it has a sizeable uncertainty of the order of 40 x 10~ n |2(J (or maybe 
even larger, as discussed in [plH). 



The latest SM predictions and the present world average of the experimental result (1) 
differ by ~ 3cr (a theoretical and experimental error combined in quadrature), if preliminary 



results of the evaluation of the vpl contribution from |13, Q are used together with the 



estimation for the lbl given in PTJ. The significant progress in the reduction of the error on 



the experimental side and the stabilization of the SM prediction of a^ l ad makes it essential for 
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the possible implications of the (g — 2) M results to be reanalysed. The issue of whether this is 
a signal of supersymmetry is already being addressed ||25|| . Here we focus on the implication 
of the latest (g — 2)^ results for the light-Higgs boson scenarios in the non-supersymmetric, 
CP-conserving Two-Higgs-Doublet-Model (2HDM) in its version called "Model II"; this is a 



continuation of earlier studies [24, 26, \2 



In sec. 2 we collect some results based on the recent calculations of af, M . As a reference 
for the vpl contribution we take the e + e~ data-driven analysis done by Jegerlehner ]T3| , FJ02, 
where the new CMD-2 results were used. The difference between the experimental data and 
the SM prediction, Aa^ = a e * p — a^ M , can be used to derive stringent constraints on the 
parameters of models, which give additional contribution(s) to a M . We calculate the interval 
Sa^, which can then be used to constrain any such contribution, at 95% CL. In sec. 3 we 
introduce the 2HDM(II), which we wish to constrain using this interval, and discuss briefly the 
existing constraints. In sec. 4 we derive the 95% CL limits on the parameters of this model, 
using <5a M . Section 5 contains the combined constraints, while the conclusions and outlook are 
given in sec. 6. 



The g — 2 for the muon: the new experimental 
and theoretical results 



Here we collect the SM contributions (and their uncertainties), which we take into account 
in our analysis. First we discuss the hadronic contributions (see table 1). We use the higher 
order contribution from and, for the lbl scattering contribution, we take an estimate from 
[^0|1 . The leading vacuum polarization contribution is taken from a preliminary result of an 
analysis by Jegerlehner, FJ02 []I3|j , where the experimental input is based only on the e + e~ data, 
including the latest ones from CMD-2 p2| . We sum all the hadronic contributions, adding in 



quadrature the corresponding errors. This leads us to the result for given in the last row 
of table 1 (which we label by the author of the analysis of the vpl contribution). 

TABLE 1 



Hadronic contribution 


in 10~ n ] 


ho 161 


- 100 (6) 


lbl [20| 


80 (40) 


vpl @ 


6889 (58) 


had [FJ02] 


6869 (71) 



We take the QED and EW terms from || and ]7|, §, respectively (see table 2). We then calculate 
the total SM prediction presented in this table, by adding the QED and EW to the full hadronic 
contributions, and by adding in quadrature the corresponding errors. This leads us to the SM 
prediction (we label it, as above, by the author of the analysis of the vpl contribution): 



[FJ02] 



116 591 726.7 (70.9) x 10 



-li 



(2) 
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Taking the new world average we calculate the quantity Aa M , defined as the difference 
between the central values of the experimental and theoretical (SM) predictions for a M . The 
error for this quantity we estimate by adding in quadrature the corresponding experimental 
and theoretical errors, o = ^Ja^ xp + o 2 SM . Next we calculate the regions of 5a^ allowed at 95% 
CL, assuming Gaussian errors. This leads to an interval symmetric around Aa M , quoted in the 
last row of table 2. 

TABLE 2 



SM contribution 


in 10~ n ] 


QED 

had[FJ02] 
EW 


116 584 705.7 (2.9) 
6 869.0 (70.7) 
152.0 (4.0) 


tot 


116 591 726.7 (70.9) 


Aa M (a) 


303.3 (106.9) 


lim(95%) 


93.8 < 5a M < 512.8 



The 95% CL interval 5a^ thus obtained is positive, and hence it leads to an allowed positive 
contribution (an allowed band). At the same time it also leads to the exclusion of any negative 
contribution to the a M . Note that the additional positive contribution to a M can be up to a few 
times larger than the EW one. 



Using results for the vpl contribution obtained by the HMNT group ||14j| , from their "exclu- 
sive" analysis, the values we get for a^ M , Aa M (a) and Sa^, are numerically very close to those 
obtained above, using the FJ02 analysis. Keeping all the other contributions as before, but 
with the HMNT(ex) results for the vpl term fll4| , we get, in units of 10~ n , 



[HMNT(ex)] : Aa fl (ff) = 297.0 (107.2) 87.2 < <5a M < 507.4, (3) 
while their "inclusive" analysis leads to a more stringent constraint, namely 

[HMNT (in)] : Aa M (a) = 357.2 (106.4) 148.7 < 5a^ t < 565.7. (4) 

Note that the above 95% CL intervals are positive for both the HMNT results, just as in 
the FJ02 case. Also note that the relative differences between the upper bounds in all three 
analyses are small (below 10%). However, the use of HMNT "inclusive" analysis leads to a 
lower bound, which is relatively much higher (up to 70%) than in the FJ02 and HMNT(ex) 
cases. 



The other recently published SM predictions [K| [12|, to which the BNL paper [Ij refers, 



were obtained from analyses of the vpl contribution based on the older, often not very precise 
e + e~ data. To improve the accuracy of the estimation of the vpl part, r decay data were 
included in addition in those analyses. However, predictions obtained there for a^ M and for 
usm are not very different from these new preliminary results used by us. It is worth mentioning 



that the following trend is observed (see discussion e.g. in [fL3"l ): if one uses the r-decay data 
in the calculation of vpl, then its estimation for the value of vpl increases while that for the 
uncertainty decreases. On the other hand, the preliminary FJ02 and HMNT analyses of the 
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vpl contribution, the ones used here, rely solely on the low-energy data for the e + e collisions. 
The accuracy of these analyses increases significantly with respect to the earlier analyses of 



the same kind thanks to an inclusion of new, high-precision measurements |22| and the use of 
more refined theoretical methods. Moreover, the central value for the vpl contribution did not 
increase in these new analyses. Thus, the preliminary analyses by FJ02 and HMNT lead to a 
deviation from the experimental value for the muon (g — 2) larger than the published results 
mentioned above, which differ by "only" 1.6 to 2.6 a, as pointed out in |J. 

Let us conclude this general part by the following comment. It becomes clear that already 
at present there is a need of a coherent and a comprehensive error analysis for all components 
contributing to the calculation of a^ M . It will be even more nessesary in a near future, when the 
BNL experiment will reach the planned accuracy. Since an estimate combining all the different 
contributions does not exist for a^ M , and hence for Aa M , which is needed for estimating the 
effects of new physics, is not available at the moment, a simplified error analysis, as presented 
in this section, is unavoidable. This is sufficient for a rough estimation of new effects; however, 
to reach a final conclusion, better error analyses are necessary. 



3 2HDM and existing constraints 
3.1 The model 

The non-supersymmetric, CP-conserving 2HDM ("Model II") based on two doublets of 
complex, scalar-fields 0i,02- This is a simple extension of the SM, in which only the Higgs 
sector is enlarged. To avoid possible large effects due to the flavour- changing neutral currents 
(FCNC), the 2HDM potential can be chosen in a ^-symmetric form, i.e. without (</>i,02) 
mixing. In the general case, the potential can have terms, characterized by a mass parameter 
/i, which break the Z2 symmetry softly; see e.g. |30], |3T |. 



The 2HDM has five Higgs particles: two neutral Higgs scalars h and H, one neutral pseu- 
doscalar A, and a pair of charged Higgses H ± . Their masses are free parameters of the model. 
Other parameters are: the angle a, which describes the mixing in the neutral Higgs-scalar sec- 
tor; tan/? the ratio of two vacuum expectation values of scalar doublets, tan/3 = Vijv\\ and the 
parameter fi. Small values of [i seem to be more natural from the point of view of the FCNC 
effects [[Kj. It is worth noticing that for such a case, a non-decoupling of the heavy Higgs sector 



can be realized |30, 31 . 



In the 2HDM one can choose the Yukawa couplings in a few different ways. Here we consider 
the Model (II) implementation, where one doublet of fundamental scalar fields couples to the 
li-type quarks, and the other to the <i-type quarks and charged leptons. This way FCNC 
processes are avoided at the tree level |29|, |28|| . This Higgs sector is identical to the one in the 



MSSM; however, in the 2HDM (II) considered by us, there are no tree-level relations between 
parameters as in the MSSM case. Therefore even for very heavy supersymmetric particles, the 
2HDM (II) and MSSM have very different phenomenology. 



5 



To be more specific, let us consider the ratios of the direct coupling constants of the Higgs 
boson h or H to the massive gauge bosons V = W or Z, as well as to the fermions (i.e. 
Yukawa couplings) for the w-type quarks and d-type quarks and the charged leptons, to the 



corresponding couplings for the SM. They are determined in terms of angles a and j3 p8| , |30f . 
For Xi = 9i/{9i)sM (and similarly for H) we have, in a form suitable for a simultaneous 
discussion of h and H: 

Xy = sin(/3 - a), Xy = cos(/3 - a), Xy = 0, (5) 
Xu = Xv + cot (3xv , Xu = Xv - cot Pxv , xi = -*75 cot (5, (6) 
Xd = Xy-tan/5xy, xf = Xy + tan /?x y , Xd = -*7s tan /3. (7) 

Here we have (Xy) 2 + (Xy) 2 = 1- Observe a pattern relation between these couplings (for h or 
H): (Xu - Xv)(xv ~Xd) = 1 - Xy, or (xu + Xd)Xv = 1 + X«Xd, found in |0| . 



For Xv — 1 a U couplings of /i have the SM values, the couplings of H to gauge bosons 
are equal to zero, while one of the couplings of H to fermions may differ considerably from 
the corresponding SM one, for a small or large tan/3. If Xy — 1 then the i7-boson has SM 
couplings, while h has very different properties: Xy = and the Yukawa coupling Xd can be 
large, for large values of tan/3. This is a case that may correspond to a light-scalar scenario 
discussed below. 

The Yukawa coupling Xd, relevant to a Higgs boson coupling to a muon, plays a basic role 
in the calculation of the 2HDM contribution to a M . It is equal (up to a factor —^75) to tan/3 for 
a pseudoscalar and H + . If in addition Xv = s ^ n (P — a) = 0, then the same holds for a scalar h; 
more precisely then |x^| = tan/3 (see eq. (7)). In the calculation of the two-loop contribution 
to a M , the coupling of H + to a scalar h is involved as well, and it is given by 



. h _ I , M h \ h , M h ~ ^ 



with the same normalization as that for an elementary charged scalar particle in the SM. For 

Xv — one gets Xh+ = (Xd ~ VXd)(^ — f j2 )/i'^^H+)- ^ e see that this coupling depends on 
the parameter /1. In this paper we consider only the case with /x = 0. A more general case will 
be studied elsewhere. 



3.2 Existing constraints 



Many searches for a light Higgs particle in the 2HDM (II) were performed at various energies 
and machines; the most systematic studies were performed at LEP. All existing LEP data, see 
for instance [27|, |32|, |33|, |34], |35| , allows for the existence of one light neutral Higgs boson, h or 



A, with a mass even below 20 GeV. According to the results presented in the left panel of Fig. 
1, the other Higgs particle (A or h, respectively) should be heavy enough to avoid the exclusion 
region in the (Mh, Ma) plane, given roughly by Mh + Ma > 90 GeV. 

This is in contrast to the case of the SM Higgs boson, which should be heavier than 114.4 
GeV (95% CL). Also the MSSM Higgs particles should be heavier than ~ 90 GeV @. An 
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Figure 1: Left: The (M^Ma) exclusion plot from OPAL [3^|. Right: The upper limit on the 
(Xv) 2 (DELPHI preliminary results [|4] 



analysis of the Bjorken process leads to an upper limit on the coupling of h to the gauge boson 
Xy- This limit, obtained at 95% CL, is presented in the right panel of Fig. 1. We see that this 



coupling is much smaller than 1 for ^ 50 GeV. The Yukawa couplings Xd of a very light 
scalar or of a very light pseudoscalar, with mass below 10 GeV, are constrained in the form 
of upper limits by the low-energy data |36|, j37H , whereas LEP experiments |35| do that for for 
masses ^ 4 GeV (see Figs. 4, 5). It is only the analysis of the decay Z — > h/A r y at LEP [27 
that gives both the upper and lower limits for \xd\- Note that \xd\ is equal to tan/3 for A and, 
if Xv = 0' a ^ so f° r h- 

The constraints from the T — > h(A)~j process, mentioned above, have been measured by a 
few groups [P^j . We present their results in Fig. 4 (lines denoted by K, N and L). Unfortunately, 
the corresponding predictions have large experimental and theoretical uncertainties, the latter 
ones both being due to the QCD and relativistic corrections. Nevertheless, as we will see 
below, the constraints coming from this process, even when accounting for some additional 
uncertainties, play an important role in closing a low mass window for the scalar h. 

Finally, note that in the 2HDM there is an important lower limit on the mass of H + , coming 
from the NLO analysis of the b -»• 37 data, given by M H + > 500 GeV at 95% CL p9|. 



4 Constraining 2HDM(II) by g — 2 for the muon data 



We apply the Sa^, obtained in sec. 2, to constrain parameters of the 2HDM (II) (see also 
earlier papers PBl lO, ED]]). We assume that the lightest Higgs boson, h or A, dominates the 



full 2HDM (II) contribution, i.e. we have a 2HDM 



"fJt ! 



or gc (a simple approach in 
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This approach should hold for masses below 50 GeV, according to results presented in the left 
panel of Fig. 1. For higher masses, which are also considered here, this is essentially equiva- 
lent to assuming a large mass gap between the lightest one, h or A, and the remaining Higgs 
bosons, which lead to a light-/i or a light-A scenario. The relevant one- and two-loop diagrams, 
studied in || |26|. 40 and [j4~T , [26|j , respectively, are shown in Fig. 2 for the h and A contributions. 





f(W+,H+) 



Figure 2: One- and two-loop (W + and H + loops are only for a /i-exchange) diagrams. 



According to the LEP limits, discussed in sec. 3.2, we assume that h does not couple to 
W/Z, and therefore we neglect the PV-loop in the light-/?, scenario. We include, however, an 
H+-\oop with M H + equal to 400, 800 GeV (and for fj, = 0). 

In Fig. 3 we present the contributions to a M obtained for an h (solid lines) and for a A 
(dashed lines), assuming Yukawa couplings Xd equal to 1. For both h and A, the one-loop 
[|, ^6|, f|0| and two-loop [§4], ^6f results are shown separately. For the purpose of comparison, 
a one-loop H + contribution is presented. The one-loop diagram gives a positive contribution 
to a M for a scalar, whereas it is negative for a pseudoscalar, independently of the value of 
the Higgs-boson mass. The signs of the two-loop contributions are reversed, these diagrams 
contribute negatively (positively) for a h (A) case [^TJ . These two-loop diagrams can give large 



contributions, since they allow us to avoid one small Yukawa coupling with a muon in favor 



of a coupling with other, potentially heavy, particles circulating in the loop |^2|, [11]]. Indeed, 
the contributions of two-loop diagrams dominate over the corresponding one-loop ones when 
the mass of h or A is above a few GeV, as shown in Fig. 3. As a result, in the two-loop 
analysis, based on a sum of the one- and two-loop (fermionic and bosonic) contributions, a 
positive (negative) contribution can be ascribed to a scalar h with mass below (above) 5 GeV 
or a pseudoscalar A with mass above (below) 3 GeV. 



In our calculation, and contain contributions proportional either to x\-, equivalently to 
tan 2 p for A and for h (if Xy = 0), or to XdXu — — 1- Assuming = 5a^ for a light-/;, scenario, 
and = 5a^ for a light-A one, and using the estimate of the interval Sa^ from table 2, we can 
derive constraints on tan (3, for h and A. They are, as expected, in the form of allowed regions 
(the area between thick lines in Fig. 4) for masses below 5 GeV for h, shown in the left panel 
of Fig. 4, and for masses above 3 GeV for A, the right panel of Fig. 4 (see also 



S 




Figure 3: The (absolute value of) individual contributions to a M from a scalar h (solid line), 
a pseudoscalar A (dashed line) and a charged Higgs boson H + (dotted line). The one-loop 
contribution for A and H + , and the two-loop one for h, are negative. Two-loop diagram 
contributions for A and h (denoted "1") are based on the down-type fermion loops only. For 
h also, results with an additional contribution due to the if + -loop are shown: line "2" ("3") 
corresponds to M H +=800 (400) GeV. 



5 Combined 95% CL constraints 

The 95% CL constraints from the (g—2)^ are presented in Fig. 4 (area between thick lines in the 
left and right panels) together with current upper limits from LEP, from the Yukawa processes 
@] (see ALEPH, DELPHI and OPAL results). Also the lower limits from the Z -> h(A)"f [27] 
can be seen in Fig. 4. In addition, the upper 90% CL limits from the T decay (lines denoted 
K,N and L, with the K results rescaled by a factor 2, as discussed in |26||) and from the Tevatron 



dc , are presented as well. 



We see that our two-loop analysis, based on the latest (g — 2) M data and on the FJ02 
estimation of a h ^ d (vpl), if combined with constraints from other experiments, allows for the 
existence, in the 2HDM (II), of a pseudoscalar with mass between ~ 25 GeV and 70 GeV, and 
tan p above ~ 25. The mass region for A allowed by (g — 2) M data, between ~ 3 and 25 GeV, is 
excluded by the constraints from LEP, based on OPAL and DELPHI data. On the other hand 
the constraints from the Tevatron close the pseudoscalar-mass window above 70 GeV. For a 
light scalar the combined constraints are even more severe; if the old constraints from the T 
decay data are taken into account, the area allowed by the latest (g — 2)^ data disappears. Note 
that the exclusion of a light h is in agreement with the conclusion of the theoretical analysis 



reported in [43 
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Exclusion 95%C.L. for h in 2HDM(II) 



Exclusion 95%C.L for A in 2HDM(II) 




Figure 4: Current 95% CL constraints for h (left) and for A (right). The (g — 2) M data 
(this analysis) give the allowed area shown between thick dashed lines. Also the upper limits 
from LEP experiments (ALEPH, DELPHI, OPAL), Tevatron, as well as from the T decay into 
h{A)^f and the rj decay (for h only) are presented. The lower limits from LEP measurements 
on Z — > /i(A)7 are shown as well. See text for details. 
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6 Conclusions and outlook 



Using the latest precise measurement of the (g — 2) M and comparing it with the improved, 
theoretical estimations of the SM contribution, we derive the 95% CL 5a^ interval to be used 
to constrain an additional contribution to a^, beyond the SM ones. It allows one to constrain 
strongly the additional contribution, which arises in a CP-conserving, non-supersymmetric 
2HDM (II) for a small parameter /i, studied in this paper. The additional contribution, allowed 
at 95% CL, has to have a positive sign, and can lead to a clear prediction for a light-scalar and 
a light-pseudoscalar scenario in the model considered here. These two scenarios correspond to 
the case when one of the Higgs bosons, h or A, is very light, much lighter than the other Higgs 
particles of the model. It should be further noted that both of these scenarios are in agreement 
with existing data from other experiments. An exchange of such light particle dominates in 
the one- and two-loop contributions to the a M . Constraints from (g — 2)^ are such that they 
exclude a light h with a mass above 5 GeV and a light A if its mass is below 3 GeV, as the 
corresponding contributions are negative in these regions. 

Our two-loop analysis presented in this paper is based on the newest (g — 2) M data and on 
the (preliminary) FJ02 estimation of a^ ad . Combining the constraints from the (g — 2)^ data 
with those from other experiments, a pseudoscalar with mass between ~ 25 GeV and 70 GeV 
and 25 ^$ tan/3 ^ 115 is allowed. However, a light scalar is excluded. 

The main results will hold also if the g — 2 constraints are based on the HMNT results for the 
vpl. For the HMNT(in) case the window for a pseudoscalar will be even smaller: 35 < Ma < 70 
GeV and tan/3 between 40 and 120. The results will not change significantly, if the uncertainty 
for the lbl contribution is up two times larger than in the estimation we used in the analysis. 
If such an error is still added in quadrature, then the lower bounds go down, with respect to 
results given in Fig. 4, by a factor ~ 1.4. There will be no visible changes in the upper bounds. 

Finally, we stress a need for a coherent and comprehensive error analysis for the SM con- 
tributions to the a M (see a very recent paper [44]). 
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